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Starting from the quantum kinetic field theory [E. Calzetta and B. L. Hu, Phys. Rev. D37, 
2878 (1988)] constructed from the closed-time-path (CTP), two-particle- irreducible (2PI) effective 
action we show how to compute from first principles the shear and bulk viscosity functions in the 
hydrodynamic-thermodynamic regime. For a real scalar field with A$ self-interaction we need to 
f^ ' include 4 loop graphs in the equation of motion. This work provides a microscopic field-theoretical 

(«_^ , basis to the "effective kinetic theory" proposed by Jeon and Yaffe [S. Jeon and L. G. Yaffe, Phys. 

^^ ' Rev. D53, 5799 (1996)], while our result for the bulk viscosity reproduces their expression derived 

^^ I from linear response theory and the imaginary-time formalism of thermal field theory. Though un- 

avoidably involved in calculations of this sort, we feel that the approach using fundamental quantum 
kinetic field theory is conceptually clearer and methodically simpler than the effective kinetic theory 
approach, as the success of the latter requires clever rendition of diagrammatic resummations which 
is neither straightforward nor failsafe. Moreover, the method based on the CTP-2PI effective action 
illustrated here for a scalar field can be formulated entirely in terms of functional integral quanti- 
_i , zation, which makes it an appealing method for a first-principles calculation of transport functions 

►^ ' of a thermal non-abelian gauge theory, e.g., QCD quark-gluon plasma produced from heavy ion 
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I. INTRODUCTION 

In a recent series of papers Il]-0], Jeon and YafFe (JY) have derived expressions for the transport functions for a real, 
self-interacting scalar field in flat space, from first principles, using the Kubo formulas ^ and the imaginary-time 
formalism of thermal field theory. A necessary step in evaluating transport functions using the Kubo formulas is 
taking the zero-frequency limit of the time- Fourier-transformed spatial correlator of the energy momentum tensor. In 
order to avoid infrared divergences which arise in taking the w — > limit, a complicated set of resummations of ladder 
diagrams had to be performed Q . However, it was observed in Q that the same integral equations for the transport 
functions obtained from the Kubo formulas can be derived from an effective kinetic theory which takes into account the 
one-loop finite-temperature corrections to the mass of quasiparticles in the scalar theory, and to the effective vertices 
for quasiparticle scattering. In their treatment, the effective kinetic theory is presented as a physically well-motivated, 
but phenomenological, theory, whose justification is taken to be the fact that it gives transport properties which agree 
with the leading-order nonperturbative results computed using thermal field theory. 

We demonstrate in this work that, in fact, the effective kinetic equations for A^"* theory proposed by Jeon and 
Yaffe are derivable from a kinetic theory of quantum fields constructed earlier by two of us |5| following the work 
of Kadanoff and Baym et al and continued by many others H . In particular their effective kinetic equations are 
derivable in our approach from the closed-time-path (CTP), two-particle-irreducible (2PI) effective action g truncated 
at four loops for the X^'^ theory. Even though the calculation of higher loop effects is necessarily involved technically, 
we feel the quantum kinetic field theory approach is conceptually clearer and methodically simpler than Jeon and 
Yaffe's approach, as the success of the latter requires clever rendition of diagrammatic resummations which is neither 
straightforward nor fail-safe. For example, the usual form of Boltzmann equation derived in, say, M assumes 2-2 
particle scattering, which conserves particle number, but bulk viscosity arises from particle nonconserving processes 
(this was emphasized by JY). When we take into account 2-4 or 4-2 processes, we need to go to 4 loop diagrams in 
the Boltzmann and the gap equations which will no longer assume the simpler form familiar in the usual derivations 
based on 2-2 processes. This extra effort is expected for tackling higher order effects but within the same context of 
the same fundamental (not effective) kinetic theory. We show that the derivation of the effective kinetic equations 
from the CTP-2PI effective action requires only the basic assumptions of kinetic theory: namely, the existence of a 
separation of macroscopic and microscopic time scales [^ and the existence of well-defined (perhaps modified) 'on-shell' 
asymptotic particle states S. The assumption of weak coupling is also necessary in order to justify the neglect of 
yet higher-order scattering terms which arise in the collision integral for the derivation of a generalized Boltzmann 
equation. 

The apparent reported 'failure' of existing kinetic equations in the literature to account for the bulk viscosity of 
the scalar quantum field is due to the fact that previous work along this line entails perturbative expansions to 
insufficient accuracy, and not to any flaw or defect of the theory. When the calculation to a sufficient accuracy is 
performed, as done here, fundamental quantum kinetic field theory produces the "effective kinetic theory" of Jeon 
and Yaffe, in particular the result for the bulk viscosity reproduces JY's expression derived from linear response 
theory. More importantly perhaps, our method based on the CTP-2PI effective action illustrated here for a scalar 
field can be formulated entirely in terms of functional integral quantization, which makes it an appealing method 
for a first-principles calculation of transport functions of a thermal non-abelian gauge theory, e.g., QCD quark-gluon 
plasma produced from heavy ion collisions. First- principles approach with a clear bearing on fundamental physics, 
involving nonlinear and nonperturbative effects, such as those employed here for this task, and elsewhere for related 
problems, are, in our opinion, essential for the successful establishment of a viable and useful quantum field theory of 
non-equilibrium processes. (For a sampling of current activities, see, e.g., M.) 

A. Statement of the problem 

Let us begin by stating in simple terms what are the transport coefficients to be computed |I(]|-[l4 . For a real scalar 
field, there are no other conserved quantities than energy and momentum. Since there is no fundamental concept of 
particle number (equivalently, since the field describes both particles and antiparticles), chemical potential must be 



^This means that the short wavelenght quantum fluctuations of the field are in a state of near local thermal equilibrium, whose 
properties vary slowly on the scale of the correlation length for the field. We also assume a nearly Gaussian initial state, so 
that low order correlation functions are suflicient to capture the dominant physical processes. 



set to zero identically in the grand canonical (equilibrium) density matrix. Equilibrium states are parametrized by 
the vector /?'' = Pu^, where u^ is a timelike unit vector {u^u^j, = — 1) and (3 = \/T is the inverse temperature. 

In the semiclassical limit, we define the energy momentum tensor T^'' as the expectation value of the corresponding 
Heisenberg operator. T''^ may be decomposed with respect to u^ as 



T'"' = pu^ 



pF'-; 



pM^ = yy^'' + u^u'' 



F^-'u^ = 



(1) 



where 77^" is the Minkowsky metric diag (-1,1,1,1) (for simplicity we shall work in flat spacetime, although a generally 
covariant formulation is readily available) , thus defining the energy density p and pressure p. Observe that there is 
no heat conduction, namely, that in the rest frame there is no energy fiux. In this case this shall be true even in 
noncquilibrium states, since there are no currents other than four - velocity to break the isotropy of space in the rest 
frame. In other words, we are forced to adopt the Landau - Lifshitz formulation over Eckart's. (See footnote below.) 
The functional dependences p = p {T) and p = p (T) define the equation of state in parametric form. The macro- 
scopic description is completed by giving a concrete expression for the entropy flux S'^, which must be consistent 
with the first law of thermodynamics dS'^ = —P^dT'^'^. Moreover, in equilibrium, vanishing entropy production 
S^ = must follow from energy - momentum conservation T^" — 0. If wc introduce the thermodynamic potential 
<j)M = 5"^ _l_ p^j^ti-v^ consistency requires $^ = pP^^ , and the identity 



dp 
Energy - momentum conservation implies the identities (recall that u^u^ = 0) 



p,t -ip + p) ur, = 0; -ip + p) u':^ + P^'^p., = 
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-u^^X, 



Since p and p become space dependent only through their temperature dependence, we may write p_f 
similarly for p. Using the identity Eq. (g), Eqs. (|^) simplify to 
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p^TT.t, and 
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where c^ = \/p,t/ P,t is the speed of sound. 

Let us now consider a near equilibrium state, meaning that the properties of the actual state remain close to that 
of a conveniently chosen fiducial equilibrium state. There is some arbitrariness in this choice, g but, following the 
so-called Landau - Lifshitz prescription |1Q|, we shall choose the fiducial equilibrium state as that having the same 
four - velocity and energy density as the actual state. More precisely, we define the four velocity as the only timelike 
unit eigenvector of the actual energy momentum tensor T'^" (assumed to satisfy suitable energy conditions), and then 
the energy density is defined as p = T^'^u^u^. Knowing p, we may compute the temperature T and pressure p of 
the equilibrium state, and thus the departure of the actual T^'^ from its value Tq'^ in the fiducial equilibrium state. 
Observe that if we write T^'^ — T^^ + ST^'^, then by construction Ui.ST'^'^ — 0. 6T^'^ is usually parametrized in terms 
of the bulk t and shear r^" stresses, as STt"" = tP^'^ + rt"" , t/^ = 0. 

Remaining within the so-called first order formalism |12,n3| , we may write the entropy production in the noncqui- 
librium evolution as 5'^ = —Pfi.i^ST^'^ (we refer the reader to the literature for a thorough discussion of the hypothesis 
involved in this formula). Decompose 
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^There are two common choices of fiducial states. In the Eckart prescription, one chooses an equilibrium state with the same 
particle current and energy density as the actual state, and reads out the equilibrium pressure from the equation of state. Thus 
in the Eckart frame there may be energy flux relative to the particle flux, which is interpreted as heat. In the Landau-Lifshitz 
prescription, the fiducial state has the same energy density, energy flux and particle number density as the actual state. Thus 
heat is read out of the particle current. 



where 
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H^-" = ^P'^^P-'- [UA.. - u^,x] (7) 
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The condition that entropy is created rather than destroyed leads us to parametrize 

r'^- = _,yi/A.-; r = -CuPp- ,/,C>0 (8) 

where ry and C are the shear and fcuZft viscosity coefficients, respectively. We wish to compute these coefficients as 
functions of the temperature and other parameters in the theory. 

B. Transport coefficients from Icinetic equations 

Since thermodynamics alone cannot determine the dependence of the viscosity coefficients on temperature, to 
proceed we must place the macroscopic description within a more basic and com pre hensive framework, i.e., kinetic 
theory, where there is a well known method to extract the transport functions ||lj,|l5| . 

The framework is a system described by a 1 - particle distribution function /. There is a known prescription to 
compute the energy momentum tensor T^" from the moments of the distribution function. In equilibrium, / depends 
only on the inverse temperature 4- vector field /3^. The starting point is the transport equation for / 
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f = Icol[f] (9) 



where Icoi is the collision integral and M^ is the mass of the particle (with possible position dependence) . / is assumed 
to be of the form f = fo + Sf, where /o is the local equilibrium distribution 

/o = ^^\ (10) 

e|/3>"| _ 1 

where /3[^ ~ m^/Tq and 5f is a perturbation. Since the collision integral vanishes identically for local thermal 
equilibrium, we can write the collision integral as a linear integral operator Q acting on Sf 

Icoi[fo + Sf] = Q[Sf] (11) 

On the other hand, if we neglect Sf on the left hand side of the transport equation, we can write it as some differential 
operator acting on /3^ , thus obtaining a linear integral equation for Sf 

Q[Sf] = QE[dt,d,]{f]l) (12) 

The Q operator satisfies four constraints, which follow from energy - momentum conservation, namely 

^e{p°)S{no)p''Q^o (13) 

(27r) 
where flo = p'^ + Af ^ enforces the on-shell condition. Thus the equation for Sf requires four integrability conditions 



{2.t 



e{p'^)s{n^)p^^QE[dud,]{l3l) = o (14) 



The integrability conditions reduce to a system of differential equations for /3° , which are in fact the conservation 
laws for the energy - momentum tensor T^^'^ eqs. (ph. These equations allow us to eliminate time derivatives from the 
transport equation, which simplifies to 

Q[5f]^QEm{Pl) (15) 

On solving this equation, we determine the correction Sf to the distribution function, and thereby the correction to 
the energy momentum tensor. In general the terms containing Sf will contribute a term Sp to the energy density; 
thus we define the physical temperature T from the condition p(T) — p(Tq) + Sp, or equivalently Tq = T + ST, 
where ST = —Sp/p^x- Knowing the temperature T we may compute the energy momentum tensor T^'^ in the fiducial 
state, subtract it from the physical T'^'^ to determine the nonequilibrium part JT'"', and then read out the viscosity 
coefficients by matching it to the form given in Eq. M). 



C. Transport functions in quantum kinetic field theory 

From the discussion above, we may identify the main steps involved in computing the transport functions, namely: 

1) Find a description of the system in terms of a 1-particle distribution function /, and the corresponding transport 
equation, 

2) Find the structure of equilibrium states, including the expression of conserved currents in terms of /, and the 
equilibrium equation of state; 

3) Solve the linearized transport equation to obtain the response of the system to gradients in the hydrodynamical 
variables, and read out the nonequilibrium stresses. 

Step 1 is done in detail in pi, where the self - energy is computed to 2-loops accuracy, giving as a result that, for 
a A0* type interaction, the transport equation for / is simply the relativistic Boltzmann equation for Bose particles, 
with the only modification of allowing for a variable mass as in Eq. (g) . This Vlasov type correction takes into account 
the fact that the physical mass M^ of the particles is connected to the temperature through the gap equation, and 
thereby A-P becomes position dependent if T is. The (only) conserved current T'^" is defined as the expectation value 
of the corresponding Heisenberg operator, and the hydrodynamic variables are read out from it, so step (2) does not 
present great difficulty. 

The problem arises in step (3), because the Boltzmann collision operator satisfies, besides the four conservation 
laws associated to energy - momentum, a fifth constraint 



'^'^P nf„0 



{2nf 



p°) S (r!o) Qn = (16) 



associated to the conservation of particle number in binary collisions. There is therefore a fifth integrability condition, 
and the system of macroscopic equations for Tq and w'^ becomes overdeterminated. One could hope that the fifth 
constraint would be true just as a consequence of the other four, but we shall show below that in an interacting theory 
this is not the case. 

Continuing on this route, the linearized transport equation built out of the Boltzmann collision operator is not 
integrable, and the calculation grinds to a halt. If we are going to compute the bulk viscosity out of quantum kinetic 
theory, then the collision operator cannot be just Boltzmann's derived from 2-2 collision processes. A generalized 
collision operator including particle number changing terms besides the usual binary scattering terms is needed. Thus 
the fifth constraint has to be lifted to eliminate the inconsistency. However, in the "effective kinetic theory" of Jeon 
and Yaffe these new terms are not derived but rather postulated to match an independent calculation of cross sections 
from linear response theory. We feel that it is conceptually and methodically more gratifying if these terms can 
be derived ab initio from a kinetic theory of quantum fields. This is indeed possible, as our present work aims to 
demonstrate. 

Since these particle number changing interactions are higher order in the coupling constant (for pure Xcj)'^ theor 
they appear at A^th order), it is to be expected that they may be retrieved by simply carrying the calculation in 
to a higher loop order. However, there appears a matter of principle: if we are going to work to high (eventually, 
arbitrarily high) order in perturbation theory, we cannot assume that the Green functions will look anything like those 
of the free theory. Thus we must first confront the need to provide a non perturbative definition of the 1-particle 
distribution function, (which should of course reduce to the one used in [p| at first order in perturbation theory). 
In equilibrium, this problem is solved by the Kubo - Martin - Schwinger (KMS) theorem ( |16[), which implies the 
proportionality of the Fourier transforms of the Hadamard and Jordan propagators (see below). Off equilibrium, 
following IgI, we shall define the 1-particle distribution function from the ratio of the partial Fourier transforms of 
these propagators. Familiarity with the KMS theorem and the Kadanoff - Baym equations should not blind us to 
the highly nontrivial nature of this definition. With this we will then have come to a full circle of deliberations for 
consistency. 

For the specific goal laid out for this investigation, the main technical difficulty lies in the analysis of the collision 
term giving rise to the bulk viscosity, as it is due to the particle-changing processes which even in the leading order 
already involve 4 loop self-energy diagrams. This is one of the main tasks we need to overcome. 

D. Summary of the paper 

The outline above provides us with a step by step guide to computing transport functions in quantuin kinetic 
field theory, which we shall execute in the following sections. As noted above, the first step is the precise definition 
of the 1-particle distribution function, which is discussed in Sec. 2. In Sec. 3 we derive the transport equation. 



For simplicity, after showing that to lowest non trivial order the Boltzmann collision operator is recovered, we shall 
write down only the terms related to particle number changing interactions. Section 4 is dedicated to studying the 
equilibrium states of the field, with the aim of finding the precise equation of state. The results of Sections 3 and 4 
amount to a first-principles derivation of Jeon and Yaffe's effective kinetic theory from quantum field theory. Finally, 
in Section 5 we go through the actual calculation of the bulk viscosity, which in the appropriate limit reproduces JY's 
estimates from linear response theory. 

II. NONPERTURBATIVE QUANTUM KINETIC THEORY 

Our specific goal is to show that by consistently extending the existing methods of quantum kinetic field theory 
(see e.g., g) to four- loop order, it is possible to account both for the shear and bulk viscosity of an interacting scalar 
field, as computed by Jeon and Yaffe. We shall consider a purely quartic interaction, although for the application to 
gauge theories cubic plus quartic would seem closer to what is needed. Since bulk viscosity entails particle number 
changing scattering, and these processes appear for the first time at O (A'') , we must push the calculation through 
to five loops in the closed time path (CTP) two - particle irreducible (2PI) effective action (EA) [|l^,Q, which will 
yield four loops in the equations of motion for the propagators. We shall assume that the background field vanishes 
identically, so we shall look at the 2PI-EA as a functional of the propagators alone |jl9;pQ| . 

A. The model 

Let us begin with the classical action for a quartically self-interacting scalar field in Minkowski space. Using a 
modification of DeWitt's notation in which capital letters denote both spacetime {x^) and time branch (1, 2) indices 
pll , the action can be written 



S = ^cI>^DabcI>'' + Si 



(17) 



Dab = [ZbO ~ ml] cab; 



Si[(t)] = -J XbCABCD4>^(, 



.BaC±D 



(18) 



where m^ is the bare "mass" of the field, Af, is the bare coupling constant and (/)^ is the scalar field. With the benefit 
of hindsight, we shall put the wave function renormalization factor Zi, — 1, but it should be generally included. The 
two- and four-index objects cab and cabcd are defined by their contraction into the scalar field. 



CABcl^^ij'' 



= d 



[0^1 



<AV^] 



(x) 



(19) 



CABCD'/'^ '/''''/''''/''' = / rf^X 



T - m" (-) 



(20) 



We shall use cab and its inverse c to raise and lower indices, and with the use of the Einstein convention of 
summing over repeated indices, their appearance may be implicit. 

We wish to derive an effective kinetic description of this theory valid at arbitrary temperature, for sufficiently weak 
coupling A, in the case of unbroken symmetry. This assumes that the expectation value of the Heisenberg field operator 
<&// vanishes. The two-point function < <^h {x) <&// (y) > is the lowest-order nonvanishing correlation function for the 
space of initial conditions with which we are presently concerned. Therefore, let us couple an external, c-number, 
nonlocal source Kab to the scalar field as follows. 



S[4>]~ 

and construct a quantum generating functional 



5M + ^Kab^^'P'', 



(21) 



Z[K] 



D(j)exp 



Aj,B 



5M + -KAB^b""^ 



(22) 



whose functional power series expansion contains all the n-point functions of the theory. The generating functional 
of normalized expectation values is given by 

W[K] = -ih\n Z[K]. (23) 

Now, we define 

oKab 
and construct a new functional F [G] which is the Legendre transform of W [K] , 

r[G] = W[K]^^KABG^''. (25) 

It follows immediately from the above definition that 

Jg^--2^^b- (26) 

and r obeys the integro-differential equation 



T[G] ^^ifiXn 



11^,,^ 1 (5r[G] A:B t.r^AB' 



i^^exp { - ( 5[0] - -^(^^^^ ^ nC^^) 



(27) 



By expanding F in a functional power series in Ti, this equation can be solved p8[ . The solution has the form 

r[G] = -yTrlnG+ yi^AsG-^^ +F2[G], (28) 

where the functional T2 is —iTi times the sum of all two-particle-irreducible diagrams with lines given by TiG and 
vertices given by the quartic interaction. 

The functional r[G] is the two-particle-irreducible (2PI) effective action whose variation with respect to G gives the 
equation of motion for the two-point function. Because we are interested in computing transport properties of this 
theory, we will need to include those terms in the perturbative expansion for r[G] which will contribute to the bulk 
and shear viscosity in the weak-coupling, near-equilibrium limit. In order to have binary scattering of quasiparticles 
in the effective kinetic theory, we shall need to have a term with four propagators, which appears at O(A^). In order 
to have number-changing processes such as two quasiparticles scattering into four (and vice versa) , we need to include 
a term with eight (six asymptotic on-shell propagators and two internal lines B), which appears at O(A^) in the 2PI 
effective action. This means taking into account the Feynman graphs in Figures 1 to 5 123]. Taking the functional 
derivative with respect to Gab yields a formal equation for the two-point function of the theory, 

hDAB-in{G-^)^g-]^TAB+l^AB = ^ (29) 

where we have singled out the tadpole term T, Fig. 6. The remainder of the self energy (which we shall refer to 
as the self energy, for short) is given by the sum of the graphs in Figures 7 to 11 (Note: observe that in the graphs 
Figs. 1-4 all internal lines are equivalent; in Fig. 5 we have instead two sets of equivalent lines, marked a and h in 
the Figure. Thus this last graph gives rise to two different graphs upon variation, i.e.. Figs. 10 and 11). This is just 
the Dyson equation for the inverse propagator, where the self energy is already expressed in terms of the propagators 
themselves. In the sense of ||l^,Q we say the self energy has been slaved to the propagators. 

There are two ways we can proceed. We can either right-multiply G^^ ^ or left-multiply G'~^^ into the equation, 
obtaining the right-multiplied and left-multiplied Dyson equations, respectively. Note that only the tadpole term in 
the is invariant under simultaneous translations of the A and B spacetime indices. The higher-order terms all violate 
translation invariance in the equation of motion for the two-point function as a consequence of slaving - they describe 
the dissipative processes by which the system approaches equilibrium [g. From now on we shall set S = 1. 



B. Nonperturbative properties of the propagators 

Our strategy is as follows. In equilibrium, the propagators are translation invariant, and their Fourier transform 
are simply proportional (KMS theorem ||lq| ). Out of equilibrium, we write 

G(x,x')= /-^e^P"G(X,p) (30) 

J (27rj 

with u = X — x'and X — {x + x') /2. We assume that G{X,p) is slowly varying with respect to the center of mass 
variable X. 

Before we start, it is useful to display the properties of the propagators which actually follow from their definition 
as path ordered products of field operators. We shall consider 8 different propagators, namely, 

a) the four basic propagators, appearing in equation (p4): Feynman G^^ =< T {<^h {x) <^h {x')) >, where T stands 
for time ordering, Dyson G^^ =< T (<&// {x) $h [x')) >, where T stands for anti time ordering, positive frequency 
G^^ =< ^H {x) ^H {x') > and negative frequency G^^ =< ^h {x') ^h (x) > .The Feynman and Dyson propagators 
are even. We also have 

G" = G22*; g12 = G^i*; Gi2 {x,x') = G^i {x',x) (31) 

As a consequence, G^^ and G^^ i^iP) are even functions of momentum, while G^^ {X,p) = G^^ {X, — p) .Moreover, 
G^^and G^^ {X,p)a,re real, and G^^ iX,p)* = G^^ ■ Finally we have the identity 

G" + G^^ = G^2 + G^^ (32) 

which follows from the path ordering constraints 

G" =eit~ t') G^i +9it' - t) G^" (33) 

G^^ = e{t- t') G^^ + 6{t' - t) G^^ (34) 

b) The Hadamard propagator Gi = G^^ + G"^"^ =< {^h (x) , ^h {x')} > is real and even and therefore also is Gi {X,p). 
The Jordan propagator G = G^^ — G^^ =< [^h (x) , ^h {x')] > is imaginary and odd, and so G {X,p) is odd but real. 

c) The advanced and retarded propagators 

Gadv [x, x') = -iG (x, x) 9 {t' - t) , Gret (x, x' ) = Gadv [x , x) = iG (x, x) 9 {t - t') (35) 

or else 

Gret = * (G" - Gl2) ; Gadv = ^ (G^^ - G^^) (36) 

Once Gret is known, we can reconstruct G as 

G (x, x) = (-i) [Gret (a;, x) - Gret [x , x)] (37) 

So 

G {p) - {-l) [Gret [p] - Gret {-p)] = 2Im G,et (p) (38) 

where we have used that Gret {x,x')is real, so Gret {—p) — Gret (p)* • Also observe that Gadv [p] — Gret ip)* 
Since the retarded propagator is causal, it satisfies the equation 

Gret — 9 {t ~ t) Gret (39) 

And therefore the real and imaginary parts of its transform are Hilbert transforms of each other 

Gret {P) = ^ I '^ ^ ■ Gret (^,p) = ^G.et {p) + ^PV f -^Gret (^,p) (40) 

ZTT J p^' — U) + le 2 ZTT J p" — U) 



ReGret (p) = -PV / ^ImG.et {lo,p) (41) 

■K J UJ — p^ 

This implies in particular that the real and imaginary parts are orthogonal to each other 

I diO ImGret {cJ,p) ReCret (w, p) = (42) 

All other propagators can be decomposed in a similar way. For example, since 

G" {X, X') = I [Gi (X, X') + G (X, X') sign (t - t')] = ^ [Gi {x, x') - i {Gret {x, x') + Gret {x', x))] (43) 



G" {X,p) = ^ [Gi - 2zReG,et] ; G^^ {X,p) = ^ [Gi + 2zRcG,e*] (44) 

To give a nonperturbative definition of the one particle distribution function /, which is the focus of attention in 
quantum kinetic theory, we shall assume that the partial Fourier transforms of the Hadamard and Jordan propagators 
are proportional 

Gi == sign (/)[! + 2/] G (45) 

Introducing a density of states function A (p) 

G(p) = 2^sign(/)A(p) (46) 

then 

Gi = 27r [1 + 2/] A (47) 

G^i = 27r [0 (/) + /] A = 2nF'^^A (48) 

G^^ = 27r [e (-/) + /] A = 27ri^i2^ (49) 

In equilibrium, / is the Bose- Einstein distribution function (KMS theorem) . It can be assumed Eq. ( |4q ) serves as 
the definition of the function /, valid to all orders in perturbation theory. Observe that, since the relevant Fourier 
transforms are distributions (e.g., in free theories), this definition may only be applied if both Fourier transforms have 
the same singularity structure, which in last analysis is a restriction on allowed quantum states. In what follows, we 
shall assume these restrictions are met. 

C. The nonperturbative retarded and Jordan propagators 

In the approximation where only terms linear in the gradients of the Fourier transforms of propagators are retained, 
it is possible to write down a nonperturbative (in the coupling constant) expression for the retarded and Jordan 
propagators. 

Let us obtain an equation for Gret from, say, the equations for G^^and G^^, namely 

Gret = I (G" ~ Gi2) ^ ^ (G^i - g22) (50) 

leading to 

— 1 = DGret ~ -zT-llGret + ^relGreti (51) 

where we have used that T12 = 0, and 

n,et = nn + ni2. (52) 



Next we perform the Fourier transform. Since we are only interested in computing the transport coefficients, we only 
need to keep terms which are first order in gradients ( this approximation which is formally invoked in the derivation 
of kinetic theory may not be always useful when dealing with realistic physical conditions, see, e.g., [^). Therefore, 
in computing the transforms, we shall drop all second derivative terms. The free term D = O — m'^ transforms into 



D — —p + ip 



d 1 



□.Y - mi 



dX 4 
We drop the D'Alcmbcrtian as they contain second derivatives: 

D p2 + ip— - ml 

The tadpole (for a generic propagator G) reads in position space 

TiBG-AG"(x,a;)G(x,2;') 

We write 

A 



-Gi (x.x) G [x.x') 
4 '^ ' ^ ^ ' ^ 4./ (2^f (27r) 



' ^ '\#^e-G,(x,,)G(X,p) 



and retain only terms linear in gradients to obtain 

^Gi(a;,x)G(x,x')-T / ^^e'^"G (X,p) 



u 

G^iX,q) + - — G^{X,q) 



The contribution to the equation has the form 



-SM'iX)-- 



t d {6AP) d 



2 dX dp 



G{X,p) 



where 



Let us write the remaining term collectively as 

UG = Jd^yU{x,y)G{y,x') 
which transforms into 






{2^Y {2nf 
Keeping only first terms in gradients, this transforms into 



G(^,P 



dS 



i4 i4 

{2nf {2n)' 



U{X,q)G{X,p) + ^L_^{dxU){q)Gip) + ^U{q)dxG{p) 



and then into 



'^y f JLJLJ^lJL Jp{y-^')e^9ix-y) 



d'^y 



{27:r {2^y 



n {X, q) G {X,p) + '- idxH) {q) dpG ip) - '-d.U (?) dxG ip) 



(53) 



(54) 



(55) 



(56) 



(57) 



(58) 



(59) 



(60) 



(61) 



(62) 



(63) 
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which contributes a term 



- {{dxHret) dp - {dpilret) dx) + ^ret {p) 



Kjrr 



to the equation of motion. 

Introducing the Poisson bracket 



we may write the equation for Gret as 



{/, 9} = dpfdxg - dxfdpg 



-l--nGret + ^{n,Gret} 



where 

n=p^ + M^ - Uret {p) , M^ ^ml + SM^ 

and we get the exact (formal) sohition 



Grpt — 



Zb {P + isf + m2 - n,et (p) 



-1 _ 1 



7mpO-+0+ 



where 



{p + ief = -{p° + ief +f 



(we have displaced p° into the complex plane to account for the retarded boundary conditions). Now write 



Then 



Gret - ReG.et + ^G 



_ -2ImVL _ -sign (/) ImQ. 

G — -7^ — ; A — ■ 



\m 



TT\n\ 



(64) 

(65) 
(66) 
(67) 

(68) 
(69) 

(70) 
(71) 



D. Equation for the negative frequency propagator 



The equation for the negative frequency propagator is 



Recall that 



DG^^ _ 1tibG^2 _^ nisG^^ ^ 



^22 _ ^12 ,-r< ^* _ 1 _ " 



After Fourier transforming, we obtain 

= -O 



Lt 7T- 



^l^'G^^l + ^i^^'Hi.} 



(72) 



(73) 



(74) 



In keeping with the stipulation to consider only first order corrections to local thermal equilibrium, we shall neglect 
all terms containing both derivatives and radiative corrections. So the equation is equivalent to 



= -r^ 



G 



12 »ni2 



ir!| 



^{^,G-} 



(75) 



To separate this equation into real and imaginary parts, we must notice that the combination i\li2 is actually real 
(see Appendix). 
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E. The unperturbed theory 

The unperturbed theory is obtamed by neglecting the O (A^) terms in our equations. The unperturbed equations 
are 



and 



nQ^p^ + AP 



A{p) =S{p^ + Af 2) + O (a2) 



(76) 



(77) 



These are exact solutions of the above equations. 

Concerning the distribution function, the real part of Eq. (75) above shows that G^^ is concentrated on the zeroes 
of Q, as required by eq. (45), and, since {Qq, A} = 0,thc imaginary part becomes the unperturbed transport equation 



= A 



^^'l'-^''^ 



which is in the form of a Vlasov equation. 



III. THE TRANSPORT EQUATION 

The nonperturbative (in the coupling constant) equation we have derived for G^^, plus the decomposition eq. (|4£ 
lead in a straightforward way to the transport equation. Neglecting {51, A} as before, we write Eq. ([75|) as 



o = -n 



G 



12 «ni2 



m 



TA{Ren,F^^} 



(78) 



Since ImO = — |0| G/2, its imaginary part reduces to 

0== A 



Fi^Inir! - ^ni2sign (p°) + i {Ren,F'^} 



which is the Boltzmann equation. To simplify it even further, we observe that since sign (p°) = F^^ — F^^ , 

Imn = -IrnRret = -ImHii + iU^ = (i/2)(ni2 - n2i) 
(see Appendix), so 



= A 



i{Ref7,Fi2}-^(ni2F2i-n2iFi 



(79) 



(80) 



(81) 



This equation is formally valid to all orders in the coupling constant. However, it is convenient to consider the loop 
expansion of 11 to reduce this equation to a more familiar form. 



A. The colhsion term 



In this subsection we shall consider the expansion of the self energy 11 in terms of Feynman graphs of increasing 
loop order, as a means to obtain a definite expression for the collision term in the kinetic equation (p^). Since we 
have the relationship 1121 (p) = 11x2 (— p) (see Appendix) it is enough to analyze only the expansion of ni2. Physically 
this means considering only the gain processes, which produce a particle within a given phase space cell. The collision 
term is then obtained by subtracting the loss processes, which remove a particle therein. 

The first term in the expansion is the single two loop graph Fig. 7. To this order. 



ni2 {x, y) = -^ A^Si^ (^, y)^^\' [G^^ {x, y)] ' 



(82) 
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In momentum space, dealing with the propagators as if they were translation invariant, 

J (^TTJ (^TI") (^""j 

and, using the definition cq. (Q), we get 

^^ ^ ' J (27r)^ (2^)^ (27r)^ ^ ) \ ) K > y ) y J 

If we substitute A by its lowest order value Aq — S (jp' + Af^), this yields the collision term given earlier in [^. 
This represents binary collisions, which conserves particle number. For the reasons discussed in the Introduction, it 
leads to an inconsistency when one tries to compute the bulk viscosity coefficient. 

The first correction to En. (pj), within the two- loop theory, comes from the radiative corrections to the density of 
states, as given by eqs. ( [7l| ) and (|80|). We write 

where Sj^ is the lowest order result just discussed, and ( writing fJos = s^ + M^ for short) 

.l2.„^ _o..o.^3 /- rf Vsign (r") d^^^ (^p^) ^4^^ (^^ ^) 
(27r)^ (27r)^ (27r)^ 

Gret (r) Gadv (r) Fi2 (r) F^^ (s) F^^ (i) Pai - U12] (r) 



JEi^ (p) = 3z (2^)^ / ^'3 3^' ^ '^^^^3"-^ "'^"^"3"*^ ^ (p _ ^ _ s - (86) 



We use again eq. (^ to get 

5S12 ( ) ^ -A^ (27r)^ /• dS<5 (r>o,) dH6 {not) d^u6 (r>o J d\,5 (»o.) d^u;^ (»o^) 
2 J (2TTf (27r)^ (27r)^ (27r)^ (27r)^ 

5(p + M + w + w + s + i)cr(M + u + u;) 

^21 (s) ^21 {t) [Fl2 (u) Fl2 (^;) ^12 (^„) _ F^^ (ti) F^l („) ^^21 (^)] 

where 

a (r) = sign (r") Gret (r) G„<i, (r) F^^ (r) 

A successful contribution to the gain part of the collision term describing scattering of two into four particles (this 
being the simplest particle number non conserving process in this theory) must involve, besides the factor (1 + /p) 
already explicit in eq. (|8l|), five other factors /or 1 + / evaluated on on-shell momenta adding up to p. But 
equation ( pTj ) cannot contain a term like this, because of the interference between the two terms in brackets. After 
all cancellations, we are left with radiative corrections to the already known binary collision term. We conclude that 
to order A^ there are no contributions to a particle number nonconserving collision term arising from the setting sun 
graph. We must consider instead the higher loop graphs, figs. 8 to 11. 

B. Higher loops 

Generally speaking, we expect the collision term to describe both particle number conserving (2 — > 2) and changing 
(2 -^ 4) scattering. Because of parity, we do not expect transitions between an even and odd number of particles. 
The 2 — > 2 scattering is already present in the two-loop theory, and any further correction to it will not contribute 
to the transport functions. So from the three and four loop contributions we shall seek only terms related to 2 ^ 4 
scattering. 

Since we only seek the lowest order contribution to the bulk viscosity, we may substitute the density of states A by 
a delta function concentrated on mass shell, so the notions of on and off shell recover their usual meaning. It is then 
possible to ascertain from the momentum flow in the graph whether the condition of five on-shell momenta adding 
to p may be fulfilled: this is just the question of whether it is possible to cut the graph by going across five internal 
lines p4| . The three loop contribution cannot satisfy this criterium, and we shall not analyze it further (it only 
renormalizes the binary scattering amplitude). For the same reason, we discard the graph in Fig 9, and concentrate 
on the graphs in Figs. 10 and 11, which a priori pass the test. 
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The complete contribution to II12 from the graph in Fig. 10 reads 

— iA* 12 f d^q d^r d^s d'^t d'^u d'^v d^w 



■i J (27r)'' (27r)^ (27r)^ (277)"* {2nf {2nf {2tt)'^ 

d{q + r + s — p)S{t + u + v — q)d{u + t + r — w) 
{2G" (g) Gi2 (r) G'^ (s) G^^ (i) G^^ (u) G^^ {v) G^^ {w) 
-G" {q) G" (r) Gi2 (s) G" (i) G" (u) G^^ («) G^^ (w) 
-Gi2 [q) Gi2 (r) Gi2 (s) g22 (i) g22 (u) G^^ („) G^^ (^i,)} 

In a true contribution to 2 — > 4 scattering, the six on - sheU momenta involved (including p) are irreducible, in the 
sense that there are no other linear relations among them than overall momentum conservation. If we look at the 
three terms in curly brackets in eq. (pq), we see that in the second term the three momenta s, v and w are on-shell, 
but they satisfy the linear relation s + v + w — p = 0, irrespective of the other momenta. Thus this term is not 
irreducible, and does not contribute to 2 — > 4 scattering; it is another radiative correction to the binary collision term. 
The same analysis disposes of the third term, since here the on-shell momenta q, r and s are constrained to satisfy 
q + r + s— p = 0. We will disregard these two terms. 

The graph in Fig. 11 contributes 

—iX'^ 12 f d^q d' 



(2^) 



4 ' "' J {2T:f {2t:)^ {2T:f {2TTf {2T:f {2T:f {2TTf 
5 {q + r + s — p) 5 [t + u + V — q) 5 [t + r + s — w) 
{G" (g) Gi2 (r) G^^ (s) G^^ {t) G^^ {u) G^^ (v) G^^ (w) 
+G12 (q) G" (r) G" (s) G^i (t) G22 (u) G^^ (v) G^^ (w) 
-G" (g) G" (r) G" (s) G" (i) G^^ (u) G^^ («) G^' (w) 
-Gi2 (g) Gi2 (r) Gi2 (s) G22 (i) G22 (u) G22 («) G22 (w) } (89) 

Only the first term in curly brackets is irreducible. Retaining only the irreducible contributions from both graphs, we 
get the prospective particle number nonconserving collision term as 

—iX'^ ,„ -.4 f d*r d*s dH d'^u d'^v _ , , 

— — 27r) / J J J J jdir + s + t + u + v-p) 

4. ' ' J (27r)^ (27r)^ (27r)^ (27r)* (2^)^ 

a^ (-P, r, s, t, u, v) Gi2 (r) G^^ (s) G^^ (t) G^^ (w) G^^ (v) 

a^ =G^^{p + r + s) [2G22 {p + v + s) + G^^ (p + u + 1;)] (90) 

It is clear that only the totally symmetric (as a function of r, s, t, u and v) part cr^ of a^ contributes to the integral, 
so we shall assume that a^ has been symmetrized. 

To reduce eq. (M) to a more familiar form, let us assume that p^ > 0, and restrict the integral to future oriented 
momenta (that is, when a momentum is past oriented, we reverse its sign). Because of momentum conservation, 
some momenta must be future oriented, but because they are all on mass-shell, they cannot be all future oriented at 
the same time; the number of future oriented momenta can only be 4, 3 or 2. The terms with three future oriented 
momenta describe 3 — *■ 3 scattering, which conserves particle number, so they are not related to the bulk viscosity. 
With these considerations we finally get the particle number nonconserving collision term as 

rd^r0{r^)s{n)d^se{s^)sin)dH0{t°)s{n)d^u0{u°)6{n)d^ve{v°)s{n) 

^"'^~ J J^nf (2^0^ (2^0^ (2^0^ (2^0^ 

{RlSl [(1 + /p) (1 + U) (1 + /.) (1 + /*) fufv - (1 + fu) (1 + A,) fpfrfsft] 
+R2S2 [(1 + fp) (1 + /.) fsftfufv - (1 + fs) (1 + /*) (1 + fu) (1 + fv) fpfr]} 



where 



5A^ 55,4 

Ri = —{27t) cr'^{~p,-r,-s,-t,u,v); R2 = — {2tt) cr'^ {-p,-r, s,t,u,v) (92) 

Si = S {p + r + s + t — u — v) ; 62 = S{p + r — s — t — u — v). (93) 
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IV. THERMODYNAMICS FROM QUANTUM KINETIC THEORY 

Our goal in this Section is to investigate the thermodynamic and hydrodynamic properties of a quantum field, 
particularly the equation of state and the speed of sound. Our starting point is the on - shell Boltzmann equation 
([79|). To render the Poisson bracket manageable, we keep only the unperturbed fi, equation (|76|), where M^ is given 
self consistently by Eqs (BOh and mTt) , namely 



The kinetic equation can be written as 



l + fiX,p) 



(94) 



^{noJ} = Icoi{X,p) (95) 

where / satisfies the constraint 

^e (/) S (f^o) P'-Icoi {X, p) = (96) 

(27r) 

which expresses energy - momentum conservation. Our concern is to investigate this (only) conservation law, but 
first, we need to express the gap equation eq. (M) in terms of finite quantities. 

A. The gap equation 



Let us write the gap equation as 



where 



M2 = mg + m^ + y M| (97) 



M^= I ^H^o)f{X,p) (98) 



^^ ^ "^'''.SiQo) (99) 



' 4 J (2ttY 

This second quantity is actually divergent, so to evaluate it we need to regularize it first. We shall use dimensional 
regularization, writing 



ml = —fi' / i-^ ^ ' ' (100) 

"2 J {2nf P^ + M^ -is ^ ' 



where the dimensionality d = 4 — e. We also go to euclidean momenta, p*^ — > ip^ , so 

2 \b , f d^P 1 

We obtain 



Write 



2'' J 


(^2nfp^ + M^ 


XbM^ 


f A//2 \ ""/' r [1 + §; 


167r2 


V4^-"V £[!-§! 



(101) 



£|i±| = l + i,l-.) + ..... ,102, 

15 



(7= .5772...). We get 



2 XbM 



,t2 



1 , / M2 

z m 



2 \4TTpi 



(103) 



We go back to the gap equation and write it as 



ZM^ = rnl + ^ 



In + M?- 



IBtT^ \4'K\l 



(104) 



where 



We now renormahze the bare coupHngs 



Z = \ 



zXb 
16^ 



(105) 



m^ = Zm ; At = ^A 



(the actual couphng in d dimensions being /i^A;,) to obtain the physical gap equation 

M^ = TO^ + -M? 



(106) 



(107) 



where 



Mi = 



M2 ^ / M2 



■In 



M^ 



^ 167r2 {iTTfi^J ^ 

With these we turn our attention to the energy - momentum tensor. 



(108) 



B. Energy - momentum tensor 



To define the energy - momentum tensor, we write the effective action in a general curved background, and then 
use the customary definition pq] 



rpflU 



2 ST 



(109) 



where only the derivative with respect to the metric in the first time branch is taken. The effective action itself is 
given by Eq. (pq). The first term Tr In G does not depend on the metric. Written in full, the second term reads 



d^x 



{ V^^ (ni^^ - ml) G" (x,x')[,^^ - (1 - 2)} 



(110) 



As usual 



5^/^ 1 



SgtJ^iy 



-99' 



5gpa 



= -g'^'g"" 



(111) 



and so the contribution from this term to T^^ is 



-d^'d" + -T?^"- (n, - ml) 



G^\x,x')\,__ 



(112) 



In the third term, the metric appears through the ^—g factors multiplying the coupling constants. Therefore the 
contribution to T'^^'takes the form 



^''"^s AT— ^2 = - ^^'^^ [G" {x, x)] ' + A.r;'"' 



5_ 
'(^Aiiii 



(113) 
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where A^ contains all the higher order contributions. To the accuracy desired, Ab is position independent, and we 
shall not analyze it further. Adding the two nontrivial contributions we get 



T^"^ (X) = - 



dt^d-' - Irjf^-a, 



G"(x,a;') 



2 ' 



— G" 



G" + Abv""" 



(114) 



To write the first term in terms of the distribution function, observe that dx ^ ip+ 2^x- We must neglect second 
derivative terms, and observe that terms involving pdx eventually vanish because G^^ i^jP) is even in p. So 



Tf"" (X) 



d*p 



P>" - Iv'^'p' 



G"(X,rt--7?^'^ 



.2 I ^r" 



G" + Abv'''' 



We are entitled to use the unperturbed approximation for G^^ 



G' 



+ 27rS{no)f{X,p) 



p2 ^ ]\,p _ ^£ 

The expressions that appear in T^'^are divergent and we must regularize them. Let us consider 



Ti 



:flU 



ddp [pMp^ _ i^M.p2] _/^d-2)\ ^r d 



jyr 



P 



P 



{2TTf p'^ + M^ -ie V 2d / J (^2TTf P^ + M^ - ie 

We rotate the integral into the euclidcan domain and compute the integral in rf = 4 — e dimensions, so finally 






327r2 



1 1 / M2 
z 7 — — In 



4 2 V47rM^ 



(115) 



(116) 



(117) 



(118) 



We also need 



G"(x,a;) = —SM^ ^ Z-' 
As 



, ,9 „zm 



Therefore 



4 



ml+'-^G''{x,x) = Z 



4Z I ^ 87r2 



Z 



9 A / , ^, zm' 



O(A^) 



(119) 



and 



6^11 



-^G 



G" = -m" 



87r2 



- Mn - - 



zm^\^ zm^M 



87r2 



2A/r2 



47r2 



-Mi 



O(A^) 



So far, we get 



Ti^^ + -77^ 



2 , As ^11 
^B + —G 



C = -rj"" { 



\ / ? \ 2 4 4 

A / zm \ zm m 



8 V87r2 



327r2 2A 



+ - 



2A 1287r 



M^ ,M^^O(A^) 



(120) 



Next call 






d*p 



p'^p''27TS{no)f{x,p) 



(121) 



and observe that 



d^p 
{2nf 



pf^p- _ _^M-p2 



2ttS (f7o) / (X,p) = T^" + -r^'^'M^M^ 



(122) 



so 
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rriflU 



^0 



rpflU 



rpflJy 



(123) 



where 



til' ui/ 



-'o — 'I 



327r2 



m 
2A 



(124) 



T,^" = -A/77'"'; 



Af = 



A'f^TO^ Af 



2A 



1287r2 



M 2 m2 



(125) 



Here, T^'' is infinite, but state independent and conserved. It belongs to the theory of the renorniahzation of the 
gravitational action (see |2|], |g^ and references therein), and we shall not consider it further. 

Consistency requires that we actually neglect the O (A^) terms in Ay, or at least that we consider them as a true 
(temperature independent) constant. Then we can establish the following identity, which will be useful later on. 
First write 



M| = - (M' 



2 2\ 

m ) 



AP , / AP 
m 



167r2 



Air^i^ 



A 



/ 



A 



M^ M^ AP f NP 
m 



1287r2 



2 A 647r2 yA-Kfj? 



constant 



(126) 



Then observe that 

a) Kf depends on temperature only through the physical mass M2, and 
b) 



dAP 



.f2 



A 



M 



M 



ji 



327r2 



■In 



/ AP 
\47r^2 



-1 9 



(127) 



This is the identity we need below. Thi s expression for the energy momentum tensor is equivalent to that given 
by Jeon and Yaffe. In particular, Eq.(127) implies that energy momentum conservation follows from the transport 
equation. 



C. Entropy flux and the H theorem 



Let us mention also the entropy flux 



5'' = 2 



(27r) 



e (pO) p'^27t5 (f^o) {(1 + /) In (1 + /) - / In /} 



Associated with this, entropy generation is given by 

d'^p 



Si^f. = 2 



i2nY 



/) 2TrS (r!o) 



In 



(1 + /) 



/ 



/, 



col 



The positivity of this integral expresses the H theorem. Let us write 

Pol = -^2^2 + h^i 
where the first term is the usual binary collision term 

d^re (rO) S {ft) d-^sO (5°) ,5 (17) dHO (i") 5 {n) 
h^2 — o- 



{2ttY {2ttY {2ttY 

5{p + r^S-t){{l + fp) (1 + fr) fsft - (1 + /.) (1 + /*) fpfr} 



(128) 



(129) 



(130) 



(131) 



and the second term invo lves the number changing interactions, already given in eq. (91) 
When inserted in Eq. (129), we find 
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S^^ — H2^2 + H2^i 



(132) 



where ^2^2 is the usual result [Q 



H2^2 



In 



(27r)^ (2^)^ 



(2^)' 



(2^)' 



<5 (p + r - S - t) {(1 + /p) (1 + Jr) fsft - (1 + /.) (1 + ft) fpfr} 



{l + fs){l + ft)fpfr 

whereas [from /2^4 in Eq. ( pi| 

1 f d'^pO {p") 5 {n) d-^re (r") S {ft) d*s9 {s°) 5 (Q) dHO {f) 5 {n) d'^uO {u^) 5 {n) d*v9 {v°) S {ft) 



H. 



2^4 



(27r)^ 



(2^)' 



(2^)' 



(2^)' 



(2^)' 



(2^) 



{R1-R2) 



In 



{l + .fp){l + fr){l + fs)[l + ft)fufv 



5{p + r + s + t — u — v) 



{l + fu)il + fv)fpfrfsft 
[(1 + fp) (1 + fr) (1 + fs) (1 + ft) fufv - (1 + /„) (1 + /.) fpfrfsft] 

is new. Thus the H theorem demands the inequality 

Ri>R2 



(133) 



(134) 



(135) 



We expect that the integral will be dominated by grazing collisions, where one of the reactants and one of the products 
carry essentially all the momentum. In this limit, Ri ~ 2i?2 (see cq. (^2|)), so the H theorem is satisfied. 



D. (Local) Thermal equilibrium states 

Our next concern is to investigate the equation of state, for a local equilibrium state described by a Planckian 
distribution function /q as in Eq. dlQ). The energy momentum tensor is decomposed as in Eq. (1^). The thermal 
component T^'' admits a similar decomposition 



d^p 



T^^T = / 7^ P>''2^<5 (r!o) /o {X,p) - ptu^u^+ptP^'' 



where 



{2.r 



i4 

PT^ I TTTI {np)''2T:5{no)fo{X,p) 



{2nr 

Since pT and il/^ are scalars, we may compute them in the rest frame 



Pt = ^r duj 



M 



^l3iAj _ 2 



\/uj^ - A'P 



(136) 



(137) 






duj -5 y/uj^ - Af 2 

2 n,f2 



For the thermal pressure, we find 3pT — Pt = —M^M^, so 

PT = 3 {pT - M^M^) 
The total energy density and pressure are then 

p = PT + I^f] p^PT-I^f 



(138) 



(139) 



(140) 



Th e eq uilibrium entropy fiux takes the form Sq — pjd^ — Tq'^ j3,y — {p + p) f]'^ — {px + pt) P^ ■ On the other hand, 
Eq. (|12|) yields S^ = $^^ - T^^Pu, where 
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(2^ 



'(/) p'"27r(5 (fio) In 



1-e-l^-P^ 



(141) 



This form of the thermodynamic potential brings to our attention other equivalent expressions for the thermal 
pressure 



and 



-1 /'C30 

T TT Jm 

Stt^ Jm 



(142) 



(143) 



Observe that Eqs (127) and ( |142| ) imply the thermodynamic relationship Eq. (0) (here and henceforth, we shall use 
d/dT to denote a total temperature derivative, that is with respect to the explicit temperature dependence through 
/o as well as the implicit dependence through JVP. We shall use d/dT when we mean only the former). Indeed, Eq. 
(142) implies 



rj^dpT 



M^ dA-P 



But pt + Pt ~ P + P, and 



^dp_ ^j.dpT_ ^dAf 
dT dT dT 



(144) 



(145) 



So Eq. (g) follows from (127). This concludes our study of the local equilibrium states 



V. LINEARIZED TRANSPORT EQUATION 

Under local thermal equilibrium, the transport equation is violated. We have I^oi = 0, while the transport part (for 
/ >0) 



' A'P 

dXt' 2 '>" dp^ 



/o = /o (1 + /o) 



P^p'Pu,. - -MlP'^ 



(146) 



Recalling the decomposition Eq. (0) and assuming the macroscopic equations Eq. (m , the LHS of the transport 
equation becomes 



/o (1 + /o) 



^P^.P.H'^'' - ^ \ ip.uf 



^TM% \ u\ 

3 2 'T ( ,\ 



(147) 



This plays the role of the Qe differential operator in Eq. ( [l5| ) 



A. The linearized collision term 



At this point we need to shift our attention to the right hand side of the transport equation, Eqs. (13C), (131) and 
(pi|). The collision term vanishes identically under local thermal equilibrium, so we need to consider a distribution 
function deviating from it. Write 



/-/o + /o(l + /o)x 



(148) 



Since Icoi [/o] = 0, only the deviation contributes to the collision integral. We keep only linear terms, and write, by 
analogy with Eq. (147) 



5Icol = fop (1 + fop) [Sh^2 + «2^4] 

where, upon introducing the momentum space volume element 



(149) 
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Dp 



d^pe{pO)s{n) 

— — 3 /Op U + Jop) 



(2^) 



we have 



and similarly 



(5/2^2 = CT / DrDsDt 6 {p + r - s - t) 



{-Xp - Xr + 2xJ 
[(1 + fop) (1 + /o.) /o./ot] 



(5/2^4 = / DrDsDtDuDv 



RiSi 



-Xp - 3xr + 2x«] 



' [{1 + fp) {1 + fr) [l + fs) [l + ft) fufv] 

where Ri^2 and (5i^2 were defined in Eqs. (§2, [9^). 



i?2'52 



-Xp - Xr + 4xs] 



[(1 + /p) (1 + fr) fsftfufv] 



(150) 



(151) 



(152) 



B. The method of moments 



Given Eqs (147) and (14£), the linearized transport equation can be rewritten as 



1 u\ , „ 

-p.p.H^'' - -^ { ip.u)' 



3 



M^ 



'-TM% 



K[X] 



(153) 



where K is an hermitian operator in the space of functions defined on the positive energy mass shell. We further 
introduce an inner product in this space by defining 



(^,x)= Dp^*{p)x{p); (x) = (i,x) 



(154) 



For our purposes, it will be enough to forfeit a rigorous solution, and to seek instead a solution using the method of 
moments. This entails first writing Eq. (153) in the orthogonal basis built ou t of the monomials I, p^, P^p'^, etc. 
(always with respect to the inner product Eq. (154), with Dp defined as in Eq. (150)), and then truncating it to only 
the first few moments. 

To simplify our notation, let us adopt the local rest frame, and write lo = p^ = —u.p. Let xo = 1 be the first 
element of our basis. The remaining functions are (i = 1, 2, 3) 



Xi=^ 



(1) 



X2 =^ 



(c.2) (^) - (..3) ^1) ^^3) ^^) _ ^^2y 



(l) (^2 



(1) (C.2) - (^)^ 



(155) 



91 =p; 92 =p 



{u;f) 



m 



To this we must add five independent functions built out of the binary products p^p^ (there are only five independent 
functions, because j? = uP' — M^ is not independent of the above). The simplest procedure is to think of these 
monomials as the composition of two spin 1 objects; the spin zero component of the composition is precisely pp' , and 
the spin 1 part, being antisymmetric, will vanish, so our functions are the five I — 2 spherical harmonics. For example, 
calling p± = Pa: ± ipy , we may choose 



Ym = {p\,p+Pz,pI ~p+p-,p-Pz,pl) ; 

We also have the relationships (sec Appendix) 

dp [1 - 3c2] 



2 > TO > -2 



M = T^ 



dT 



M2 - \TMjj, 



(156) 



(157) 
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,,_^,dp [M^-lTM^Tc's] 



NP 



iTM^T 



dT 



(iuf) = 3clT^ 



In terms of the new functions, Eq. (|153|) reads 



-r:;Y^w^ + r\{oj)x2 + c 



(w)Xi 



where 



[m2 - ^TM}r] 

'^^ dT 









(158) 



(159) 



(160) 



(161) 



W hen we expand the operator K [x] we notice that, when truncated to the subspace spanned by the functions eq. 
( |l55| ), the operator matrix acquires a block form, with one block corresponding to the Xa functions (a = 0, 1 or 2), 
another to the q^, and yet another to the Y^ functions. Since there are no g^ functions in the left hand side of Eq. 
(|160|), we may as well write 



X=^b^Y"-~-r[A + Bxi + Cx2] 



(162) 



Since K [lo] — 0, the B coefficient will remain undetermined ( the left hand side of Eq. ( |160[ ) is orthogonal to lo, so 
the system is integrable). We will set B = for the time being, and postpone further discussion until wc enforce the 
Landau - Lifshitz conditions. 

To determine the bm coefficients, we must solve the linear system 



T"'H'^ = X™"6„ 



where 



\^mn 



{Y"',K[Y'']) 



By symmetry, the X matrix must be diagonal 
(for the positivity of h, see Israel ^M) leading to 



6> 



tm _ _Y'>njTiJ 

5 y 



(163) 
(164) 
(165) 
(166) 



To find the A and C coefficients let us expand 



K[l] = 



{K[l]) 



1- 



A'M^^ 



(^)Xi 

;^,xi) 



(^,Xi> 



■/? 



X2 



X2 
I) 



(x 



+ 7/ 2 
(X2 



(167) 



where we have used (lj, K [x\) = 0. If only binary scattering is considered, then also {K [x]) = and {K [1]) = /3 = 0. 
In general, then, {K [1]) ^ f3 ^ j. Therefore these equations admit an approximate solution with C ~ 0, yielding 



Y — —T™HY 



Co; 



Co 






where we have used {K [1]) < 0, as follows from the inequality Eq. (135) and the identity 

DpDrDsDtDuDv {R2 - Ri) 5i 



{K[1])^AX^ 



[{l + fp){l + fr)il+fs)il + ft)fuM 



(168) 



(169) 
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C. The temperature shift and the bulk stress 



It can be seen from Eq. (168) that the correction to the distribution function has two components. The one 
associated with the H^'^ tensor contributes to shear stress, but it does not induce a change in the energy density, 
and therefore it is compatible with the Landau - Lifshitz matching conditions. The constant shift of x by cq, on the 
other hand, affects in principle both the energy density and the thermal mass Mt ■ So, to enforce the Landau - 
Lifshitz conditions, it must be partially compensated by a temperature shift. Concretely, if we call T the temperature 
of the fiducial equilibrium state, such that p (T) is equal to the energy density in the nonequilibrium state, then 
the temperature appearing in the local equilibrium distribution function /o must be To = T + 6T. The effect of this 
temperature shift is the same as that in the coefficient B in Eq. (162). 

The distribution function and temperature shifts in turn produce a shift SM^ in the physical mass, which likewise 
does not affect the transport equation. However, both ST and SAP are relevant to the bulk stress. Observe that there 
is no shift in the four velocity u'^. 

The three displacements cq, 5T and SAP are related by the constraints that the gap equation must hold, and the 
total energy density in the nonequilibrium state must be the same as in the equilibrium state. Write the gap equation 
as 



AP 



f 



(M^i 



= ^M^ 



(170) 



The linearized equation then reads 



As a matter of fact, 



So finally 



1- 


1 


XdMl' 
2 (9M2 




A 


['a^^+-H 



dT T2 



5Af 2 = M\8T + M^ Co 



(171) 



(172) 



(173) 



where 



Ml. = T'^M\ 



(174) 



Since the gap equation is enforced, we can look at the (cosmological) constant A as a function of M^ , and 

-1 



2 Jn,r2 



bkf = -^m4^sm 



(175) 



then 



Sp — p^T^T 



Actually 



dpT 1 



- -M4 



DAP 2 



dpT 1 



M^co + (uj^) Co 



^ -M4. - ^^ 
5M2 2 ^ 2T 



so 



(176) 



(177) 



Sp — p.tST + 
And since the total energy remains the same, 

P.tST = -Co 



(,,2^ _ \^TM}. 



(^2^ „ \11tm^^ 



Co- 



(178) 



(179) 
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Let us apply the same reasoning to the bulk stress, which results from both the departure of the pressure from 
p (T) and the direct contribution from the new terms in the distribution function 



T = c^ptST 



dpT 
DAP 



1 



M^ 



Mlco + \ [{lo^) - M^ (1)] CO 



Now 



so 



dAP 



-^M 



#2 



-Co 



^^"3 



{-') + 



- -±TMi 



Using Eqs. (|l6|), (H^) and (jTH), we get 



3^5^ 



l(A'[l])| 



(180) 



(181) 



(182) 



(183) 



D. Shear stress and viscosity 

The shear stress can be read directly out of the new terms in T!^'^ . In the rest frame, we get 



from which we can read out the shear viscosity 



R 

Iff 



(184) 



(185) 



To estimate 77, it is enough to keep only the leading (binary scattering) contributions, so 77 ~ A ^ . On dimensional 
grounds, h ^ T^ and R ^ T^ , so we recover the usual result, ry ^ T^/A^. 

E. Built viscosity 

As expected, things are not so simple with the bulk viscosity. We can read it out from Eq. (|183|) as 



[AP - '^TM%]' {{co^) (!) - (co^) {uj)y 



3T' 



m 



mm 



(186) 



However, in evaluating it we must consider that (1) is logarithmically divergent in the massless limit, so we must 
correct the sheer dimensional estimate to (1) ^ T^ In (M/T). As for the size of |(i4r[l])|, observe that the integral 
is dominated by the Rayleigh - Jeans tail, where /o ~ T/u; > 1. Thus \{K [1])| -- X'^T^F (AP). Since the overall 
units are Mass'^, it must be \{K [1])| ^ T^ /A'P. For the remaining elements we may use the conventional estimates 
(^3) r-^T^^ pr^ T"*, and thus obtain 



c 



AP_ 

A4T3 



n 2 



AP - -TAP 



:T 



In^ {M/T) 



(187) 



which reproduces JY's Eq. (5.6) Q. 

In the limit in which the bare mass vanishes, or equivalently in the T — > cx) limit, we may write on dimensional 
grounds 
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M^ - -TM'jj. = -fJ-M^^ - \M^ (188) 

and since M^ ^^ AT^ itself, Eq. (187) reduces to C ^ XT^ In^ (A), again in agreement with JY 0. 
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VI. APPENDIX 

A. ni2 is purely imaginary 

To show this, observe that ni2 is the sum of all IPI Fcynman graphs with two external vertices, one carrying a 1 
index and the other a 2 index (this follows from it being the result of opening one 12 leg in each 2PI vacuum bubble). 
It can also been represented as 

n - ^^' 



(189) 

b=Q 

where Fiis the usual (IPI) effective action, and (p"- the background field. The effective action has the structure 

Fi = i y dxdy m (x) D (x, y) {0} (y) + i [0] (x) N (x, y) [0] (y)} + O (0^) (190) 

where {(/)} = (p^ + cf)'^, [<f>] ^ (f)^ — (fp'] both D and N are real, N is even, and D is causal. Therefore 

^1^ = lJdy{D {x, y) {0} {y) + [0] (y) D (y, x) + 2iN (x, y) [0] (y)} + O (0^) (191) 

and 

ni2 = -iN (x, y) + ]^[D (x, y)~D (y, x)] (192) 

The real part of ni2 is odd, and its imaginary part even, which shows that its Fourier transform is purely imaginary: 
Write 

ni2 (x, y) = I -^ e'P(--^)ni2 (p) , (193) 

then the identity 11^2 (a;,y) = —1112 (y,a;) becomes indeed II^j [p) = — II12 {p) . 

We may use the same argument to find that 1121 {x,y) — ni2 (jj^x), so 1121 ip) = ni2 {—p) is also imaginary. We 
also find 

Hii = iN (x, y) + ]^[D (x, y)+D (y, x)] (194) 

from where Inillii {p) = («/2)(ni2 + 1121). 
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B. (lu) and (cj^) 



Our objective is to compute 



1 /"^^ 



Recall the identity 



w/o (1 + /o) - T 



dT 



(195) 



This and Eq. (H) may be used to establish the identity 



Similarly 



(w^) - Af 2 (w) = ST (p + p) = 3T' 



dp 



2 "P 2 



dT 



(196) 



\ / 2 '^^ '^ dT 



(197) 



So 



3X_^2rfp[^^'-i™TC?] 



(^3) ^ t2 



dT 



M 2 - iTAf ?p 



(c^) = T- 



dp [1 - 3c2] 



dT 



M2 - iTM2j, 



(198) 
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VII. FIGURE CAPTIONS 

Figl. Two loops contribution to the CTP effective action. 

Fig2. Three loops contribution to the CTP effective action. 

Fig3. Four loops contribution to the CTP effective action. 

Fig4. Five loops contribution to the CTP effective action. 

Fig5. The other five loops contribution to the CTP effective action. Observe the two sets of inequivalent lines, 
marked a and b. 

Fig6. One loop contribution to the self energy (tadpole graph) 

Fig7. Two loops contribution to the self energy (setting sun graph) 

Fig8. Three loops contribution to the self energy. 

Fig9. Four loops contribution to the self energy 

FiglO. Another four loops contribution to the self energy. Cutting as shown, we go across five internal lines. The 
symmetric cut also goes across five lines. 

Figll. The final four loops contribution to the self energy. Cutting as shown, we go across five internal lines. 
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